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' The effective gg-interaction is derived from Lagrangian QCD in the front form by 

means of the flow equations. It coincides with previous results. 



Before plunging into the paraphernalia of field theoretical details, Franz Weg- 
^ ■ ner might allow me to outline his method [1] in yet an other language. Working 

. with a Hamiltonian in an explicit representation, it is always possible to di- 

^ . vide the complete Hilbert space arbitrarily into two pieces, the P- and the 

i Q-space. Suppose, one is unable to solve the eigenvalue problem, because of 

^ ' say computer limitations. The method of flow equations allows then to uni- 

tarily transform the Hamiltonian matrix analytically into a block-diagonal 
■ effective Hamiltonian, i.e. 



[qHP QHq) ^ ^^^^ I QH,sQ 



where P and Q = 1 — P are the respective projection operators. (The meth- 
ods also works for Hamiltonians with more than 2x2 blocks.) The 'reducible' 
5^ I matrix can then be diagonalized separately for either of the two spaces. The 

transformed Hamiltonian is a function of the continuous parameter < / < oo, 

j^H{i) = m,H{i)], v{i) = mi),H{i)]- (2) 

The generator of the transformation r){l) is subject to some choice [3]. The 
Hamiltonian is separated conveniently into a block-diagonal part Ha{l) and 
into the rest Hr{l) = H{1) — Hd{l) which is purely off-diagonal: 

(PH{l)P \ z.n^( PH{l)Q\ 

^''^'''-y QH{l)Q)' ^'^''> ~ \QH{l)P J-^"^' 

In many cases of practical interest one can interpret this rest as a 'resid- 
ual interaction'. If the rest vanishes, or if it is exponentially small, one has 
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solved an important part of the problem. The generator is always off-diagonal: 
Pr]{l)P = Qri{l)Q = 0. The flow equations (2) for the diagonal and the off- 
diagonal block can then be divided into two coupled equations, 

JliPHP) = PrjQHP - PHQr]P, j^{PHQ) = Pr]QHQ - PHPr/Q, (4) 

and into the trivial identity PrjQ = PHP HQ — PHQHQ. The corresponding 
equation for QHQ will not be needed. 

A possible measure for the 'block-off-diagonality' is O = TtPHQHP > 0. Its 
derivative is determined by the flow equations, 

j^O ^Tr[PriQ(QHQHP - QHPHP)) (5) 
+ Ti({PHPHQ - PHQHQ)Qr]P) = 2Ti{Pr]Qr]P) < 0, 

because rj is anti-hermitean. Since O is positive with a negative slope, the 
generator of Eq.(2) ensures a monotonic decrease: The off-diagonal block tends 
to vanish in the limit I ^ oo. With increasing / the effective Hamiltonian ifgff 
becomes 'more and more block-diagonal'. 

Almost simultaneously with [1] and with a similar objective, Glazek and Wil- 
son [2] have proposed the 'similarity transform' to render a Hamiltonian 'more 
and more diagonal'. 

One should note the subtle but important difference between the two meth- 
ods. Glazek and Wilson address to generate a band matrix whose width A 
is decreased by analytical unitary transformations. In the asymptotic limit 
A ~ 1/Z — > the Hamiltonian becomes diagonal. It is asked for too much, per- 
haps, to find a transform which diagonalizes a Hamiltonian completely. Flow 
equations have been applied to many problems, including solid-state physics, 
as well as field and gauge field theory. I do not even attempt to review all the 
work but refer to Wegner's talk [3]. I have taken great inspiration from his 
work with P. Lenz [4] . 



1 Applying the flow equations to quantum chromo-dynamics 

In the next two sections I apply the flow equations to Non-Abelian gauge 
theory in the front form along the lines of earlier unpublished work [5]. The 
thinking in terms of block matrices is almost ideally suited for gauge theories, 
where the Fock-space expansions give quite naturally a block structure of the 
Hamiltonian, cf. for example Fig. 2 in [6]. It is thus reasonable to identify 
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the P-space with the lowest Fock-space sector, the one with one quark and 
one anti-quark {qq). The Q-space is then hterally the 'rest'. In front form, the 
Hamiltonian proper is H = P~ , with 

p- ^T + V + F + S. (6) 



This as well as all other formulas in this section are taken without further 
notice from the Compendium [7], using the same notation and conventions as 
there. The kinetic energy T is the only piece of P~ which survives the limit 
g ^ 0. The coupling constant g is related to the (strong) fine structure con- 
stant as = g'^ / {Anhc) . The vertex interaction V is the rclativistic interaction 
per se. It is linear in g and changes the particle number. The instantaneous 
interactions F and S (the 'forks' and the 'seagulls'), are quadratic in g and are 
artifacts of the light-cone gauge. Both may change the particle number. - As 
a technical trick and for gaining more transparency the instantaneous interac- 
tions are first omitted. They will be re-installed at the end of the calculation. 
By the same reason all color factors will be suppressed first and re-installed 
at the end. The Hamiltonian is then 

H^T + V. (7) 



Since the front-form V has non-vanishing matrix elements between Fock states 
only, if the particle number differs by 1 and only 1, the initial P-space Hamil- 
tonian is diagonal, PHP — T, see below. As opposed to that the effective 
Hamiltonian in the P-space has an effective interaction U, 

PH,sP = T + U. (8) 



The 'residual interaction' is then just the vertex interaction = V. The 
Q-space Hamiltonian remains a complicated operator for all /. It will not be 
addressed here. Objective of the present work is the calculation of U. 

The way the flow equations work is that V{1) de-creases while U{1) in-creases 
for growing /. As emphasized by Wegner [3], flow equations generate a whole 
series of many-body interactions. They can be classified with the power of the 
coupling constant g. For the reason of demonstrating the method, I restrict 
here and below to the lowest non-trivial order in g without further mentioning. 

The flow equations (4) are re-ordered and re-labeled as 

r,{i)^[T{i),v{i)], ^^mim], ^^m),v{i)]- (9) 
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The relativistic kinetic energy is diagonal for all I: 



+ d\dq + a\a^ . (10) 



The conventions in the Compendium include that [d^q\ stands for the integra- 
tion over the longitudinal momentum p'^ and the transversal momenta pi, as 
well as the summation over helicity, color (glue) and flavor. The single particle 
energies depend potentially on the flow parameter through the effective mass 
m^(p; /). This dependence will be suppressed below, however, since it comes in 
higher order of g. Hence forward, p+ and will be collected in the 3- vector 
p — (p'^jPl). The interaction term is 



T^/n 1 M gi f[d^ f [d%] ^(3)^^ - - 7^ 

^(0 = ^r^ ^= ^ (Pi - - Ps) 9{Pl,P2,P3]l) 



Pi ypi yP3 

b\b2a3 {uiizU2) - d\d2a-;i (?72^3^'i) + aid2h {^2^1^^) + h.c. (11) 

The three-gluon vertex is omitted here because it does not contribute to the 
effective interaction to the lowest non-trivial order. The effective coupling 
g{pi,P2,P3'J), with the initial value g{pi,P2,P3'J = 0) = 9, contains all the 
Z-dependence. The generator 77 = [T, V] becomes 



W)^-7f^ / '{Pl-P2-P3)V{Pl,P2,P3]l) 

X (616203 [M1/3M2] - 4403 [^^2/3^1] + 01^263 [v2fiU3]j - h.c. (12) 
The structure of rj is very similar to V, because T is diagonal, thus 

v(pi,P2,P3;l) ^9ipi,P2,P3;l)c{pi,P2,P3;l)- (13) 



The commutator with T generates the 'commutator function' 
c{pi,P2,P3] I) = e(pi; /) - e(p2; /) - e{p3; I), 



(14) 



with the single particle energies e{p; I) = (m^(p; /) + p±)/p'^. Symmetries like 
ViPiyPiyPi ~ Pi'J) = QipiiPiiPi ~PiJ) are so obvious that they, as well as 
^(Pi)Pi)Pi~Pij — ~^(Pi)Pi)Pi~Pij 0' ^i^^ be mentioned hence forward. 
Similarly, the derivative dV{l)/dl — [ri{l),T{l)\ becomes 



dV{l) 
dl 



Vl67r3 



[dM [[dM f[d'q3] 



S^^\pi -P2-P3) 



(15) 
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X (616203 [ui^3U2\ - 4^^203 N^3^'i] + ald2h [^2^1*3]) + h.c. 
Finally, the derivative dU{l)/dl = [^(0)^(0] becomes 



[dhi] f[d'q2] f[d\s] fW,] f[d%] f[d% 



dU{l) _ _ 1 
dl ~ ~T67r3 



Pi ypi 

S^^^ (pi - pi - Ps)^^^^ {P2 -P2- P3) b\bva3, al,dld2' [1*1^31*1'] [^^2'^'^^2] 



X 



I5 



Pl,P3]l)9{P2,P2,pi]l) +V{pLP2,P3^i)9{Pl,Pi,P3-,l)) + 



S^^\pi-Pi - P3)S^^\P2 - P2 - pi) [dld2'a3>,alb\by\ [mi^Mi'][v2'/3'^^2] 

X (^(Pl^Pl,P3;o^(P2,P2^P3^o +^(p2,P2^p^^ (i6) 



plus terms containing commutators like [a'^db, b'^d'^a] or [a^db, a^a^a] contribut- 
ing only in higher order of g. After execution of the various commutators it is 
obvious that U has the general form (c/. Eq.(8) ) 

[d'qi] J [d%] J [d'q[] J [d'q'2] 5(3) (pi - Pi + P2 - pD 

X KA{ 42^?i (l2]q'i,q2,l)- (17) 
It suffices thus to calculate the derivative of a c-number 



dU 
'dl 



^s{^(pt-Pi^) + ^(Pi^-Pt)) 



167r3 

X {f]{pi,pi,pi-pi;l)g(p2,P2,P2-P2J) 

+ 9 {pi, Pi, pi - Pi, 1)V{P2,P2,P2 -P2,l) 

uiPl, Ai)7^m(p/, a;)] d^^{q) [v{p^, \'^)Yv{j>2, A2)] 



X 



9" 



(18) 



The polarization tensor d^^ arises by the summing over the helicitics, i.e. 
d^y{q) = J2\^fi{q, ^)^l{q, The 4-momentum of the photon is denoted by 
q^. The step function is 9{x) = 1 for a; > 0, thus 9{x) + 9{—x) = 1. 



2 Integrating the flow equations and shaping the notation 



Eq.(9) gives the first order flow equations in operator form. After evaluating 
the matrix elements they reduce to two coupled algebraic equations: 



5 



j^9iPi,P2,P3;l) = -cipi,P2,P3;0) r]{pi,p2,p3;l), (19) 



?i,P2,P3;0= c{pi,p2,P3;0) g{pi,p2,P3;l). (20) 
They can be combined into a single one 

J^{PUP2,PZ] I) = -C^(Pl,P2,P3; 0) g{pi,P2:P3; I): (21) 

and integrated with the possible solution 

9{Pi,P2,P3;l) ^ 9 exp (-/ c^(pi,P2,P3;0)) . (22) 

Since the p's appear in combinations like {pl,Pi,pl — Pi, I) (see f.e. Eq.(18) ), 
I introduce conveniently the energy differences 

Wq = -c(pl,pi,p{-pi;0) = (pi-pi)~ -pi~ +pT, (23) 

Wq = -C{P2,P2,P2 - P2, O) = (p2 " P2)' - P2 + Pt ■ (24) 

Typically, they occur along the quark {wq) and the anti-quark line {wg). They 
are in simple relationship to both the 4-momentum of the exchanged photon 

?M=PiM-^^i/^ + ^MY =P2m-P2m + ^mY' (25) 

and to the Feynman 4-momentum transfers along the two lines 

Ql = -{p[ - p^f = q+wg, Q\ = -{p2- pL? = q^w^- (26) 

The latter need not be equal in a Hamiltonian approach. The null 4-vector 
r]^ is defined by 7]^ = (0, 0, 2) and cannot be confused with the operator rj. 
The Feynman-momentum transfer is always positive and a more physical 
quantity than an energy difference. The w's will be substituted therefore by 
the Q's as long as no misunderstanding can arise. In fact, I shall use often the 
mean-square momentum transfer 

Q'^\iQl + Ql)-Y^^, + ^,), (27) 

and the mean-square difference 

SQ'^l{Ql-Ql)^^-^{wg-w,). (28) 
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Also the polarization tensor d^i^ in Eq.(18) can be simplified with these defini- 
tions. Since d^y{q) = —g^u + [fJ^Qu + fly(ltj)/(l^ appears always in combinations 
with the spinors, the Dirac equation (pi — pi)fju{pi)j^u{p[) = can be used 
to write 

111 

q,u{p,, X,)ru{pl K) = -^vMpi, ^ih'uipi, X[). (29) 

In the context of Eq.(18) one can use thus d^y{q) = —g^u + V^.VvQ'^ / {<1^T ■ 

Following Glazek and Wilson [2], 1 introduce a dimensionless similarity func- 
tion /(/) (with boundary conditions /(O) = 1 and /(oo) = 0) by 

f{Pi,P2,P3;l) = ^9{Pi,P2,P3-J)- (30) 

Note that / describes the decay rate of the off-diagonal vertex interaction V. 
The rate of decay, however, is not a unique function, as will be seen below. 
In line with recent developments [3], I want to keep / = /(w; /) as a general 
function, particularly since fj can be rewritten as 

1 (d\nf{w;l)\ 

(31) 

Finally, I define 

N,^- Jdl''-^^^^f{w,; n, N,^l- N„ (32) 



which are kind of dimensionless occupation numbers. 

With these definitions I return to the problem of integrating Eq.(18). Its rel- 
evant part is 

v{pLpiJ[-Pi\'^)9{P2,pLp2-pL\^^ 

\Wg dl Wq dl ) 

The /-integration leads to 

^ dl'(r]{pl,pi,pl-pi;l')g{p2,P2,P2-P2;l')+V^^ 
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(34) 



It is natural to assume [3] that the similarity function /(w; /) is a homogeneous 
function of its arguments f{w; I) — f{w'^l), with some exponent k. Two types 
of similarity functions are considered: 



Gaussian cut-off («; = 2): f{w;l) = exp{—w^l), Ng — ^ ^ — -, (35) 
Exponential cut-off (k = 1): f{w;l) = exp{—wl), Nq — — . (36) 

Wq -\- Wq 

The first case corresponds to Eq.(22). The second shows that the effective 
interaction depends explicitly on the similarity function. The requirement of 
block diagonalizing the Hamiltonian determines the generator only up to a 
unitary transformation [3]. 



3 The effective gg-interaction from Hamiltonian flow 



Collecting terms, the effective interaction from Eq.(18) becomes 



U = g 



167r3 



Nq,N£ 



It is easy to restore the color factors. They sit only in the vertex interactions 
and go through untouched by the flow equations, thus 



jj^ 2 C^M.(g)KPl,Al)T"7^^(W,A;)]KP2^A^2)y"7'^^(p2,A2))] (Ng ^ Nq\ 
If the dni,{q) is substituted, this turns into 



167r3 



I Pi \/p'^ 



, Nq 



(37) 



The result holds actually for the general non-Abelian theory SU{nc). 
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It is equally straightforward to re-install the instantaneous interaction: 



167r3 



(38) 



cf. diagram 6*3 i in Table 5 of [7]. To the order considered here it is independent 
of /. Adding it in by U — > U + Ui gives 



= - 



167r3 



X 



IP2 

'N, Ng\ SQ^ 



(39) 



Contracting the Lorentz indices (with 7^77^, = 7^) gives 



= - 



as [ujpi, \,)T'^Yu{pI a;)] [v{pi, \',)T^^XP2, A2)] (N^ ^ N, 

Ql 



47J-2 



Pi Pt P'l P'2 



as [u{pi,Xi)T- ^Mpi, a;)] [v{p^, A'2)rV^(P2, A2)] 
47r2 



Pi pt P>\ p'2 



X 



(N, 



(40) 



One notes the non-integrable quadratic singularity {q'^ Y appearing both in the 
dynamic amplitude Eq.(37) and in the instantaneous interaction Eq.(38). In 
the perturbative calculation of the qq scattering amplitude, these singularities 
cancel each other exactly, see for example Sect. 3D of [6]. Also in the present 
approach they tend to cancel, but a residual piece remains in Eq.(40). It 
is weighted with a coefficient SQ"^ which possibly can be considered small, 
but beyond that, the term is dangerous because it violates manifestly gauge 
invariance (reflected in the matrix elements (7'''7''')). The quadratic singularity 
carries however also weighting factors which depend on the similarity function, 
see Eqs. (35,36). Inserting them, I get for the exponential cut-off: 



Nr. 



Q 



Ql + Ql' 



6Q' 



.Ql Ql 



.Ql 



Ql 



1 
0, 



(41) 
(42) 



and for the Gaussian cut-off: 
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{q^f \Ql QV 

For 5Q'^ <^ the Lorentz contracted part is obviously the same for both 
cut-offs. The potentially dangerous coefficient of (7'''7^), however, behaves 
drastically different: It vanishes strictly for the exponential cut-off. For the 
Gaussian cut-off it can be small, at least it behaves now like but its 

precise behaviour near the important region of the Coulomb singularity at 
~ must be left to future analysis. To easen my life 1 will work hence- 
forward with the exponential cut-off and its gauge- invariant result. 
Thus far, 1 have studied the Hamiltonian proper H = P~ . But it is of- 
ten advantageous to work with the Lorentz-invariant 'light-cone Hamiltonian' 
i^LC = P/iP^ [6]. Its eigenvalues have the dimension of an invariant mass- 
squared. In the intrinsic frame {P± = 0) one has i^LC = P^P^ ■ The effective 
interaction of -ffLC, is obtained simply by U =^ UP^, thus 

— ^ Cc Ai)7^m(p/, a;)] A2)7m'^(P2, Aa)] a, 1 

47r2 ^xil-x)x'{l-x') P-"' 



(Q2)2 + (5Q2)2' 
(g2)2 + (^g2)2 (^+)2 



(43) 
(44) 



where x = / P^ is the longitudinal momentum fraction of the quark. The 
color factor is Cc = Z^a^c" c'^c" C2- '^^^^ finishes my objective: the so obtained 
effective interaction for QCD coincides identically with the result of lowest 
order obtained with the method of iterated resolvents [8]. 
Acknowledgement. 1 thank Franz Wegner for many illuminating discussions. 
1 thank him also for our joint effort to shift the chaff from the wheat. 
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